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■ Abstract 

o, 

We consider the dynamics governing the evolution of a many body system constrained by an non- 
abelian local symmetry. We obtain explicit forms of the global macroscopic condition assuring that 
at the microscopic level the evolution respects the overall symmetry constraint. We demonstrate 
the constraint mechanisms for the case of SU (2) system comprising particles in fundamental, and 
Q_i' adjoint representations ('nucleons' and 'pions'). 
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^ ■ 1 Introduction and Overview 



The consideration of the influence of internal symmetries on the final state of a many body system 
begun with the pioneering work of Bethe ||]. Much of the subsequent interest in the subject arises 
from the realization that in the study of hadronic interactions and in particular in studies involving 
quark confinement, these constraints may be of decisive importance. An important progress in treating 
equilibrium systems was made employing group projection techniques. This allowed for a consistent 
treatment of abelian || and nonabelian || ||, || |(| symmetries of compact groups and a consistent 
formulation of thermodynamics of many particle systems with internal symmetries taken into account 
IH, Application of these methods to specific processes demonstrated in which circumstances the 
presence of symmetry is of physical relevance [10, 11, 12 1. 



Bowever, it is not fully understood how the symmetry-modified properties of the equilibrium 
system arise from kinetic formulation of the dynamical evolution. When an internal symmetry is 
not at work, Boltzmann's B theorem in principle assures that the statistical Bose/Fermi/Boltzmann 
distributions are the asymptotic (equilibrium) distributions, irrespective of the nature of microscopic 
interaction. Bowever, in presence of exact symmetries the equilibrium distributions are modified, see 



e.g. y, 11]. This implies that symmetry constraints introduce effective interactions of potentially far 



more complex nature than is the usual two body Boltzmann collision term. In fact it can be argued 
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that quantum symmetry constraints are the heart of the nonlocality of quantum physics. However, in 
the limit of classical Boltzmann equation evolution these are implemented by a strictly local (though 
non-linear) consideration of Fermi blocking and Bose enhancement in phase space evolution. Our 
aim in this work is to make a step towards understanding how the microscopic nonabelian symmetry 
constrains operate within the kinetic master equation description of the time evolution, leading on to 
the symmetry modified (constrained) macroscopic many particle equilibrium state. 

It is first important to convince oneself that an underlying symmetry of microscopic interactions 
does not lead in general to the desired symmetry properties of a (macroscopic) many body interacting 
system. To do this, we consider the high energy nuclear (heavy ion) collisions and specifically here 
two symmetry examples: 



a) SU(2) Isospin symmetry: The initial state transforms under a given representation of the isospin 
SU(2) group. All elementary high energy interactions are governed by the strong interaction, which 
preserves the isotopic symmetry. A final state results as a multiparticle state formed by many indi- 
vidual hadron - hadron collisions. In any of such microprocesses the isospin is conserved. However, 
proceeding 'as usual' without symmetry constrained treatment of local interactions does not assure 
that the final multiparticle state (macrostate) transforms under the same representation of the isospin 
group as the initial state, which is required for symmetry reasons. 

b) SU(3) Colour symmetry: A similar situation appears in the context of the quark-gluon inter- 
actions, especially in case that local deconfinement occurs. The initial state is a colour singlet state, 
and quark-gluon and gluon-gluon interaction, although invariant under the colour SU(3) C symmetry 
group do not assure that during its evolution a (macroscopic) many particle state, once a singlet, 
always remains a singlet colour state, which, however, it must do because of exact colour symmetry 
of strong interactions. 



As these examples show, in a dynamical (quantum) transport theory description of the approach to 
equilibrium there must exist a subsidiary condition which should be taken into account by correspond- 
ing kinetic equations governing the evolution. This condition is independent from other constraints 
related to dynamical gauged internal symmetries. For classical fluid dynamics a dynamical evolution 
equation addressing gauge symmetry has been proposed by Wong jl3|| . In many current studies of 
the dynamics of classical non-abelian fields this proposal continues till today to attract considerable 
interest [14]. However, these are constraints which have no relation to the intrinsic non- locality of the 
quantum system which we address here. 

We first note that in the case of an abelian symmetry there are no additional constraints to 
consider. Quantum number conservation on a microscopic level is fully equivalent to preservation of 
all symmetry properties on the macroscopic level. This is easily seen considering the U(l) symmetry 
related to microscopic particle-antiparticle formation: since the microscopic mechanisms produce equal 
number of particles and antiparticles (pair production), initial particle-antiparticle number difference 
is exactly preserved in the macroscopic many body state. 

Thus only presence of a nonabelian symmetry poses a true challenge. A suitable mathematical 
method how to approach this problem is identified considering the previously treated statistical equi- 
librium case. Here one decomposes a general 'macrostate' consisting of many particles into possible 
irreducible representations of the symmetry group. Then a projection technique exploiting character 
function properties of the group is used to constrain the final state. We will here use this approach 
in order to describe multiparticle evolution applying microscopic kinetic theory scheme. We will show 
that the behaviour of particle phase - space distribution functions will depend not only on properties 
of basic interaction, but that it also depends on global properties of the macroscopic system. Those 
global properties provide subsidiary constraints needed, so that the asymptotic equilibrium state has 
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properties consistent with the non-abelian constraints. 
2 The Projection Method 

Let G be a compact internal symmetry group of our system consisting of particles (objects) trans- 
forming under irreducible representations of the symmetry group. These representation are denoted 
as oti with corresponding dimensions d{pii). One denotes fj^' Vt \T,r,t) a distribution function of the 
particle which belongs to the multiplet ai of the symmetry group. Members of this multiplet are 
numbered by indexes Vi (vi = 1, . . . , d{pti)) which correspond to given values of charges related to the 
symmetry group. A subscript £ denotes other quantum numbers characterizing different multiplets of 
the same representation a. The variables (r, r) denotes a set of the phase - space variables such as 
(p, f) and t is time. 

The number of particles of the specie {a, u a , (} is: 

N £}(0 W = / dVdTf^ (r, r, t) ; (1) 

We consider a system of {Naj 1 }, (t), . . . , Na^l an (t)} particles at time t. The distribution functions 
fulfill the generalized Vlasov - Boltzmann kinetic equations, which can be written in the general form: 



%) df + v-Vf^(T u r,t) 

= E E E / dT J dT k dT l W < ^^f l) {T k ,T l ;T J ,T i ) 

aj,a k ,ai Vj,v k w Cj,Cfc,Ci 

[jffit (T 4 , r, t^'^^r, t)/g*^ (T k ,r, t)f^ {V h r, t) 

(r fej r, t)^\T h f, i)/gf i} (T h f, t)/£f i} (r, , f, t)] ; (2) 

Factors J-& v \t, r, t) are related to quantum statistics and they are equal to 1 for classical par- 
ticles, and equal to [1 ± f^ u \v, r, t)] for bosons/fermions correspondingly. Since by assumption the 
whole system transforms under given representation A of an exact symmetry group, the system under 
consideration must preserve its transformations properties during its time evolution, provided that it 
is governed by a symmetry invariant interaction. 

We now focus on the case of a quantum system and consider state vectors in particle number 
representation: , . . . ,N^2n/ ■ These vectors describe symmetry properties of our systems and 

all other variables, related to phase-space properties of the system are suppressed here. They transform 
as a direct product representation of the symmetry group G. This representation is of the form: 



ocy <S> a 2 <8> • • • <8> a n ; (3) 

A multiplicity N^' of the representation ay in this product is equal to a number of particles which 
transform under this representation: 

^ = e(e <'?< J = E N ^f ; (4) 

The representation given by Eq. (||) can be decomposed into direct sum of irreducible representations 
Afc. Corresponding states are denoted as |A^, Aa,, ; M) where AA fe is an index numbering members of 
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the representation A and J\f is a total number of particles 

* = E ; (5) 

k 

Each physical state can be decomposed into irreducible representation base states with amplitudes 
depending on phase space variables T: 

^•••^^r) = E e E ffi l A ^ Aft ;^ ; a>^ ^.(e^r); (6) 



Here appear new variables £a which are degeneracy parameters required for the full description of a 
state in the "symmetry space". 
Let us define an average weight 

This expression gives the probability that \ ■ ■ ■ , Ni^' particles transforming under the sym- 
metry group representations a%, . . . , a n combine into N particle state transforming under representa- 
tion A of the symmetry group. 

We make the statistical hypothesis that average weights (|7|) do not depend on phase - space 
variables and can be calculated alone on basis of symmetry group consideration. This also can be 
proved under the stronger assumption that in Eq. (|6|) any state with fixed A, Aa has the same weight 
(see e.g. g) 

Let us consider a projection operator P A on the subspace spanned by all states transforming under 
representation A. 

■p A 

This operator has the generic form (see e.g. |15(| ): 

V A = d(A) [ d^g)x W (g)U(g); (9) 



N™, JVftO) = E e |A, Aa; £a) (£ a ) ; (8) 



G 



Here is the character of the representation A, dfi(g) is the invariant Haar measure on the group, 
and U(g) is an operator transforming a state under consideration. We will use the matrix representa- 
tion: 

U(g)\N<Z?,...,N<%>) 
- E ■ ■ ■ ^> ■ ■ ■ D$ ■ ■ ■ D^l \ N %\ . . . , JVfc)) ; (10) 



v is a matrix elements of the group element g corresponding to the representation a. Notation 
convention in Eq. ( |T(i| ) arises since there are Ni>"^ states transforming under representation aj and 
having quantum numbers of the v aj -th member of a given multiplet. 



The statistical hypothesis identifies the average weight P ,t \ i i with a norm of the vector 



V 1 



. . . , Nv2r/- This norm can be written as 
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/iV (Ql) ••• iV (Qn) V k iV (Ql) i\r( a -)\ = |c a - a a f^All 2 - fill 



where the relation (V A ) 2 = V A was used. 

Left hand side of this equation can be calculated directly from Eqs.@ and (|io|). One gets finally 

= A^d(A) [ d»{ g )x^Xg)[D^] N ^ ■ ■ ■ pfe£]<? ; (12) 

where A'-'^ is a permutation normalization factor. For particles of the kind {a, we included in 
Eq. ( fl^ ) the permutation factor: 

^« = ^_ ; (13) 

(c) n<1 (c) ! 

The permutation factor is a product of all "partial" factors 

The permutation factor assures the normalization of state vectors: 

(N^, • • • , . . . , N^) = A {M} ; (15) 

This normalization reflects an invariance of the state vector with respect to permutations which shuffle 
indistinguishable particles. 

3 Incorporation of Symmetry 

The expression Eq. (^) is a starting point for further considerations. It provides together with Eq. ([]]) 
and Eq. @) subsidiary constraints on distribution functions f( ai ' Ui \ These conditions assure that in 
a dynamical evolution the symmetry of the system is preserved. When symmetry is conserved, then 
all weights in Eq. ( |12"D are constant in time. In a case of strong interaction and colour symmetry, all 
weights, except for the weight corresponding to the singlet state, must remain zero. 
We now convert the global constraint into a time evolution condition and consider: 

1 P {< ) ,.,<»)} = 0; (16) 
Introducing here the result of Eq. (12) one obtains: 

= ^d(A)/«W*< A ^ 

G 

+ EE -^t AW}d ^) J M9)X {A H 9 )[D^I] N ^ • • • [Dttf lZ 'og[/#;: ; (17) 



All integrals which appear in Eq. (H) and Eq. (|7|) can be expressed explicitly in an analytic form for 
any compact symmetry group. 
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To write an expression for the time derivative of the normalization factor A^^ we perform analytic 
continuation from integer to continuous values of variables Nv^ . Thus we replace all factorials by 



the T-function of corresponding arguments. We encounter here also the digamma function ip [16]: 

d logr(x) 



V>(x) 



dx 



This allows to write for Eq. (|T 



dt 



i 



E 



dt 



(18) 



(19) 



To get a consistent analytical continuation in the number of particles one should define the time deriva- 
tives dN^ a \s\/dt. We define these rates of particle number change from the integrated Boltzmann 
kinetic equation, Eq. (0), explicitly 



dN, 



(Oi) 



d,t 



E E E 

otj,a k ,oti Vj,v k ,vi QXkXl 



J dvdTjdTkdTtdTi w^v^uf l) (r k , r, ; r is r f ) 



(CO 



{ y k) {T k ^t)f^{ l \T u r,t) 



-r, 



(20) 



(CO ^'^'^(CO ! ( r ''^'*)Aci) ^^''"^fe) 

Contributions from gradient terms of Eq. (||) vanish due to Gauss law. These terms are transformed 
in surface integrals and beyond the volume occupied by the system all distribution functions are equal 
to zero. 

Eqs dl^pll 2C ) in fact constitute the global subsidiary condition which should be fulfilled by the 
microscopic kinetic equations Eq. (^). These are the necessary conditions for preserving the internal 
symmetry on the macroscopic level. Rates of change dAfjf}^/dt are related to "macrocurrents" , 
which are counterparts of "microcurrents" related directly to a symmetry on a microscopic level via 
the Noether theorem. Eq. (^) can be considered as a set of conditions on macrocurrents to provide 
consistency with the overall symmetry of the system. Therefore we believe that this equation can 
also be used as a starting point for multicomponent hydrodynamic equations with internal symmetry 
properties taken into account. One should notice that this subsidiary condition takes into account 
also surface effects for the finite volume systems. This is due to the space variables integration which 
is performed in Eqs (|l|) and (|20|). 

One easily sees that for the case of abelian symmetry the two constraints Eq. ( |T2| ) and Eq. ([17]) 
do not lead to new results: first we recall that all irreducible representations of abelian group are 
one-dimensional. Next, let basic particles have "charges" q±, . . . ,q n , and let the global charge be Q. 
Then the only consequence of Eq. ([l2] ) follows for nonvanishing weight Q = N±qi + ■ • • + N n q n , which 
is a rather obvious result. New results appear only for nonabelian symmetries. 



4 Example: Isospin 

We now consider as an example the case of the SU(2) symmetry with basic particles transforming 
under spinor (i) (fundamental) and vector (1) (adjoint) representations. This example can be realized 
by a gas mixture of nucleons and pions. To describe all group elements the three group's parameters 



a,/3, 7 are chosen in such a way that diagonal matrix elements have the well known form [15]: 
i.) for the fundamental representation (i): 



e im(a+ 7 ) cog 



m 



±- 



(21) 
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ii.) and for the adjoint representation (1): 



I>2, ± l(«>ft7) = \ ^^(l + cos/?); (22) 
45(^/3,7) = cos/3; (23) 
The Haar measure for the SU{2) group in this parametrization has the form 

27T 2n 7T 

/ d^g)f[g] = ^LjdaJd~fJ d(5 sin (3 } \g{a, (3, 7 )] ; (24) 



Any 'macrostate' is made of an arbitrary number: n n , n p , n_, no, n+; where subscripts refer to 
members of the fundamental representation neutrons, protons; and members of adjoint representations, 
7r~, 7T°, 7r + , correspondingly. Let us consider the special case when the macrostate is a SU(2) singlet. 
The weight of the singlet state is according to Eq. (^) : 



p0,0 
{n n ,n p ,ri-,n ,n + } 



2tt 2tt 7r 

JdaJdjJ dp S mPe-^ n "- n - +2n — 2 «+)(«+t) C o S ^^cos no /3 



= AW??>° v ; (25) 

where 

1Z = n n + n p + 2n_ + 2n + ; (26) 
The permutation normalization factor is here: 



.{M} _ ( n - + n Q + n + )\(n n + n p )\ _ 
~~ i i i i i > 



(27) 



n_!rio!n + !n n !n p ! 

The real nonzero values of the weight is obtained only when the argument of the exponent in Eq.( 
vanishes: 

n n -n p + 2n_ - 2n + = ; (28) 

This is equivalent to the conservation of the third component of the isospin. 

Novel behaviour is obtained only when one considers time evolution of the system. Presence of 
an exact symmetry means that the corresponding weight Eq. is constant, here we consider the 
expression: 

d-^ro - (29) 

^ {n n ,n p ,n-,no,n + } ' V > 

We note that an appropriate analytical continuation should be made. First we evaluate the integral 
appearing m Eq. (||): 



_L_ p°>° = (-i)n yv_2v ( n °\ L 

{Af} {n n ,n p ,n-,n ,n + } V > > J ft j_ l 



AIM}* {nn,n p ,n-,n ,n + } V V /_^\ "1 I ^ I ^ . ^ , ■ > (^0) 

i=0 \ / 

This discrete form is not allowing an analytic continuation which would allow for all necessary 
differentiations. However, we can write the integral also as the hypergeometric 2F1 function flGf]: 



— IrnPf' x = (-l) n °^— 2 F 1 (-n ,'R+l,n + 2;2): (31) 
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where 1Z is as defined in Eq. 



We so obtain: 



P 



o.o 



{n n ,n p ,n- ,no,n + } 

T(n_ + n + n + + l)T(n n + n p + 1) 



T(n_ + l)r(n + l)r(n+ + l)T(n n + l)r(n p + 1) T(-n Q ) f^ Q K + l + i i\ 



cosnovr ^ T(— no + 2' 



Eqs. (|25D, (|29D and (|32|), together with the condition (|28[) result in: 



op; 



,0.0 



{n n ,n p ,n_ ,no,n.|-} / dfl 



dK 



dt 



+ 2- 



dn. 



dt 



dP 



o.o 



+ 



{n n ,n p ,n^,no,n + } dno 



0n 



dt 



P 



0.0 



fjj- ~ {n n ,n p ,n-,n ,n + } ' 

The projection integrals determine the coefficients which are, explicitly: 

i0,0 



dP 



{nn,n p ,n-,no,n + } 



dn, 



-IT 



sinnovr ^2, r(— uq + i) 2 l 







E 



+ 



r(-no)^ ft+l + i i! 
cos n vr ^ r(-n + i) [^(-n ) - ^>(— no + i)] 2* 



4) no £ 



E 

i=l 



n 



K + l + i i\ 
(l + no)-^(l + no-i)] 



aP°'° , 

{n n ,n p ,n_,n ,n + } 
3ft 



fit) 



i=0 



n 



71+ 1 + i 

1 

(ft + l + i) J 



and 



dt 



+ 



+ 



dn n 

~dt 

dn_ 

~dt 
dn+ 

~dt 



dft 

(n N + 1) - ip(n n + + [ip(n N + 1) - ^{n p + 1)] 
(n* + 1) - + 1)] + ^ [VK + 1) " ^(™o + 1)] 

(n* + 1) - + 1)] ; 



(32) 



(33) 



(34) 



(35) 



(36) 



where n^r = n n + n p is the total number of nucleons and = n_ + no + n + is the total number of 
pions. Eq. ( |36| ) can be also written in the form: 



d logAW 
dt 



+ 



dn n 
dt 


riN 

E 

k=l+n n 


1 

k + 


drip 
~~dt 


njv -y 

, ? jfe 

fe=l+n p 




dn^ 
dt 


E 

fc=l+n_ 


1 

k + 


driQ 
~~dt 


^ 1 dn + 
> - H 1 

/C = l + 71{) 


fc=l+n_|_ 



(37) 



Eqs. offer the final result for the SU ( 2 ) 

case . Notably, they imply a relation for the 
number of neutral pions in a system. We thus see, that when the case of the non-abelian symmetry is 
carefully considered and not ignored, one can get relations determining also the 'neutral' members of 
multiplets. In the standard approach the multiplicity of neutral particles is obtained by introducing a 



subsidiary chemical potential which is related to the lack of chemical equilibrium of a system [17| or 
to the residual interaction with the environment 18]. 
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5 Conclusions and Outlook 

We have shown how constraints due to the preservation of the symmetry properties of a multiparticle 
macroscopic state define the path of evolution of the system. One should notice that results we 
presented are general and do not depend on the particular choice of the representation of the symmetry 
group. For different initial symmetry group representations one gets different paths, but they are all 
of a similar "shape", considering the hypothesis that the global behaviour of a macrosystem (in a 
sense of statistical physics) should not be altered if a number of particles is changed by a very small 
( "microscopic" ) amount . 

We have explicitly presented the example how our constraint works in the simplest non - trivial 
case of SU (2) symmetry group. 

Although we have studied and implemented the discreet symmetry using quantum states, whenever 
we referred here to a dynamical equation we considered the limit of incoherent state evolution described 
by the Boltzmann equation. The dynamical evolution we consider thus is described in terms of 
diagonal density matrix. This is the appropriate approach given that our main objective is to arrive 
at a dynamical derivation of symmetry deformed statistical distribution. 

We recall that quantum correlations (without symmetry) alone are responsible for the deforma- 
tion of the Boltzmann distribution into Bose/Fermi distributions, and that the Boltzmann equation 
yields this result when we allow for Fermi blocking/Bose enhancement in the collision term. In that 
line of thought, the next step would be to show that it is possible to obtain now within a dynamical 
Boltzmann equation calculation the evolution of a many body system into symmetry-deformed statis- 
tical equilibrium distribution. We are also exploring the possibility that the methods here presented 
allow the formulation of a microscopic transport theory which would obey the long range correlations 
introduced by the macroscopic quantum and symmetry constraints. 
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